G3: Genes|Genomes|Genetics

Pyruvate dehydrogenase complex (PDHc) deficiency (PDCD) is a primary mitochondrial disorder characterized by neurodevelopmental disability, altered intermediary metabolism and early mortality. Dichloroacetate (DCA), a pyruvate analogue, is a well-described PDHc activator that remains under clinical investigation for treatment of PDCD. Here, we studied the in vivo efficacy of a 5-point log concentration range of DCA on animal health and metabolism in C. elegans with feeding RNA interference (RNAi) expression knockdown of either PDHA-1 or DLD-1 homologues at graded degrees to model variable disease severity. These worm models recapitulate phenotypic features of PDCD observed in human patients, including reduced survival, delayed growth, locomotor impairment, and elevated lactate and/or pyruvate tissue levels. DCA treatment appeared well-tolerated, with no gross morphologic toxicity seen at doses up to 25 mM. Significantly improved health, survival, tissue lactate levels, and mitochondrial physiology were observed at 25 mM in pdha-1(RNAi) knockdown animals. DCA treatment in dld-1(RNAi) C. elegans models (undiluted, 1:20 dilution, and 1:100 dilution) showed significant therapeutic benefits on survival, neuromuscular function and metabolic phenotypes primarily in the moderate (1:20) and/or mild (1:100) dld-1(RNAi) deficiency strains, but not in full-dose dld-1(RNAi). Importantly, linear growth, neuromuscular activity, and mitochondrial physiology were significantly improved with DCA treatment even in the most severe dld-1(RNAi) undiluted model. Overall, preclinical modeling provides objective evidence of DCA therapeutic efficacy in C. elegans expression knockdown strains for two well-conserved homologues of PDHA1 and DLD that represent distinct genetic etiologies of PDHc deficiency, with demonstrated beneficial effects on survival, healthspan, tissue lactate, and mitochondrial physiology. These data further confirm that DCA's therapeutic effect correlates with PDHc disease phenotype severity in dld-1(RNAi) animals.

The kynurenine pathway metabolizes tryptophan into 3-hydroxykynurenine (3-HK), a precursor for ommochrome eye pigments synthesized via the cardinal (cd) gene in mosquitoes. While cd disruption was presumed neutral, we observed fitness costs in Anopheles stephensi knock-in but not knock-out cd mutants. Here we investigated this anomaly further by assessing survival, fecundity, and midgut integrity across multiple cd mutant lines. Heterozygous knock-in lines, expressing a fluorescent marker and guide RNA for CRISPR/Cas9, exhibited reduced survival post-blood feeding, larva-to-adult survival deficits, and midgut barrier dysfunction, whereas knock-outs showed no such costs. Oral supplementation with xanthurenic acid partially rescued knock-in mortality, implicating oxidative stress linked to 3-HK metabolism. Expression analyses suggest transgene insertion effects, rather than cd disruption, underlie these fitness costs. These findings highlight the importance of evaluating insertional effects in gene drive target selection and support cd as a viable target for genetic control strategies in An. stephensi.

Hearing in Lepidoptera mediates a range of ecologically important behaviours, including mate communication, predator avoidance, and acoustic signalling. In moths, the evolution of predator-prey interactions with bats has further shaped hearing through a sensory arms race, with repeated co-option of auditory organs to detect and evade echolocating predators. Despite significant prior characterization of the neurophysiology and behaviour of hearing in moths, the genetic basis of hearing is poorly understood in most insects. In this study, we identify a core set of putative auditory genes in Lepidoptera using a combination of homology-based searches from Drosophila and evolutionary rate analyses. We find 56 genes present across all species and investigate whether gene copy number varies among non-hearing and hearing lineages and among 3 different ear types. We discovered seven genes associated with ear type and one with ear presence, but did not find significant losses in gene copy number in non-hearing species. We identified three genes (btv, Dnai2, and nompB) with strong evidence of selection in hearing clades and five genes with weaker evidence of selection. We discuss the potential roles of btv, nompB, and Dnai2 in ciliary transport and the aging of hair cells, as well as the possibility of actively amplified hearing. Our study serves as a primer and resource for further gene mining and functional testing of auditory genes in moths and other insects.

Aedes aegypti is a major vector of arboviruses, and belongs to subfamily Culicinae, a diverse group of mosquitoes with homomorphic sex-determining chromosomes. Males are the heterogametic sex with a dominant male-determining locus (M locus). The M locus and its counterpart m locus are embedded in a region of suppressed recombination, with a large portion of this recombination desert showing significant molecular differentiation despite homomorphy. We developed a mathematical framework to examine M-linked genome editors that specifically target the m-chromosome during spermatogenesis, mimicking the naturally occurring sex ratio distorters (SRDs) in Culicinae that produce male-biased meiotic drives. Unlike previous models for species with heteromorphic sex chromosomes (e.g., X and Y), we incorporate features stemming from the homomorphic nature of the Ae. aegypti sex chromosomes such as varied linkage to the M locus, making the degree of super-Mendelian inheritance readily tunable. We evaluated in silico SRDs with a range of M-linkage and editing efficiencies and established the theoretical foundation for developing highly efficient SRDs that outperform several methods of population suppression. These SRDs can be tuned to mitigate impact on a neighboring population. The framework developed here is suitable for exploring SRD-mediated genetic biocontrol of pests with homomorphic sex chromosomes.

We present Nemo 2.4, an advanced forward-time individual-based simulation framework designed to model the complex eco-evolutionary dynamics and genetic basis of quantitative traits. This tool addresses current challenges in evolutionary quantitative genetics by providing unprecedented flexibility and computational efficiency. Nemo 2.4's modular architecture allows researchers to design custom life cycles by combining specialized Life Cycle Event (LCE) modules, from reproduction and dispersal to selection, crossing, and phenotype expression. The software supports diverse population models, including both Wright-Fisher (WF) and non-WF dynamics, spatially explicit models, and varying demography. Nemo 2.4 handles a wide range of genetic architectures, including both multi-allelic Quantitative Trait Loci (QTL) for general trait studies, and dense di-allelic Quantitative Trait Nucleotides (QTN) implemented with highly optimized bit-wise data structures. Crucially, it allows the simulation of QTNs on comprehensive genetic maps that incorporate other genetic elements, providing genomic-scale resolution. Key biological complexities are integrated natively: the model accommodates modular pleiotropy, dominance, and pairwise epistasis across multiple traits, facilitating the study of complex genotype-phenotype mappings. Furthermore, Nemo 2.4 models phenotypic plasticity through reaction norms and incorporates underlying liability thresholds, enabling the simulation of environmental influences on trait evolution with various forms of selection (e.g., Gaussian, linear, truncation). Due to its compiled design and memory-efficient data representations for large numbers of loci, Nemo provides a robust platform for running high-throughput simulations critical for testing theoretical predictions in polygenic adaptation and understanding evolutionary responses to changing environments.

BackgroundThe fission yeast stress-activated protein kinase (SAPK) pathway includes a conserved mitogen-activated protein (MAP) kinase cascade that regulates multiple cellular processes and is activated by several types of external stress. Understanding how Sty1, the MAP kinase in the SAPK pathway, controls these processes is complicated by the fact that different stressors can have stressor-specific effects that may be difficult to separate from the effects of Sty1 activation itself. Moreover, upon stress, Sty1 activation is usually short-lived. Previously, we developed a fission yeast strain, SISA, in which Sty1 kinase activity can be switched on in a sustained manner in the absence of external stress. This required combining multiple mutations in the SAPK pathway, including an analog-sensitive version of Sty1. When SISA cells are grown in the presence of analog-sensitive kinase inhibitors, Sty1 is inhibited, but when inhibitor is removed, Sty1 becomes hyperactive. While this strain was useful, it had several limitations. ResultsHere we describe and validate a more rationally-designed strain, SISA4, that retains the features of the original SISA strain while overcoming its limitations. SISA4 is more stable genetically than SISA, easier to use in genetic crosses, and easy to identify by phenotype or genotyping. We show that analog-sensitive kinase inhibitors 4-Amino-1-tert-butyl-3-(1-naphthylmethyl)pyrazolo[3,4-d]pyrimidine (1-NM-PP1) and 4-Amino-1-tert-butyl-3-(3-bromobenzyl)pyrazolo[3,4-d]pyrimidine (3-BrB-PP1) are equally potent for inhibiting analog-sensitive Sty1 in vivo, and we determine optimal inhibitor concentrations for converting SISA4 cells from a Sty1-inhibited state to a Sty1-hyperactive state. We also find that both 1-NM-PP1 and 3-BrB-PP1 have measurable off-target effects in wild-type cells, although these are modest and generally do not affect interpretation of experiments. Finally, using SISA4, we show that the Sty1-activated transcription factor Atf1 plays an unexpected role in maintaining cell-polarity disruption after Sty1 hyperactivation. ConclusionsSISA4 will be useful for investigating how SAPK pathway activation regulates diverse cellular processes.

Humans and other vertebrates contain two estrogen receptors (ERs), ER-alpha and ER-beta, which mediate the physiological actions of three estrogens: estrone (E1), estradiol (E2) and estriol (E3). Of these three estrogens, in vivo, E2 is the strongest transcriptional activator of ER-alpha and ER-beta, E1 is next most active, followed by E3. We studied transcriptional activation of human ER-alpha and ER-beta by E2, E1 and E3 in African green monkey kidney (COS-7) cells, which we compared with studies of estrogen stimulation of ER transcription in human em-bryonic kidney (HEK-293) cells. To our surprise, in COS-7 cells, E3 had the lowest half-maximal response (EC50) for human ER-alpha and ER-beta than either E2, which was second most active estrogen, or E1. In contrast, for human ER-alpha and ER-beta transfected into HEK-293 cells, E2 was the most active estrogen, followed by E1 and E3. Similar results were found in COS-7 cells and HEK-293 cells transfected with elephant shark ER-alpha and ER-beta. Thus, under some conditions, E3 is a more active estrogen than either E2 or E1. This suggests that E3 may be a novel physiological ligand for the ER in some mammalian cells.

Many biological characteristics arise by interactions between more than one biological organism or unit. Fertilization success in sexually reproducing species represents such an extended phenotype where both mates are required to be fertile for a successful outcome. Consequently, predictive models should account for the joint nature of reproductive performance while offering interpretable estimates for individual mate contributions. Recent advances in genomics and machine learning (ML) provide standardized, high-dimensional genetic information on one hand and computational tools capable of modeling complex biological systems on the other. Here, we construct and evaluate two-tower (TT) machine learning architectures for genomic prediction of binary reproductive outcomes and recovery of sex-specific fertility liabilities. Simulated datasets, generated under a range of genetic architectures, were utilized to compare multilayer perceptron (TT-MLP), convolutional neural network (TT-CNN), and L1-regularized linear (TT-LASSO) two-tower models. Simulation scenarios varied sex-specific heritabilities, genetic correlations, infertility prevalence, mating structure, and sex-specific infertility rates. Models were evaluated with regard to their ability to predict reproductive success at pair level and also recover true underlying genetic values for male and female fertility. Prediction accuracy increased with the underlying heritable component as expected, while sex-specific tower-scores successfully recovered latent fertility liabilities despite models being trained only on observed joint outcomes. TT-LASSO achieved the highest overall classification performance, whereas TT-MLP provided more balanced and consistent recovery of sex-specific genetic values across scenarios. An additional simulation, incorporating genotype-dependent mate compatibility demonstrated advantages of fully-connected neural networks for capturing non-additive interactions. These results indicate that two-tower frameworks provide a powerful approach for modeling reproductive traits, enabling simultaneous prediction of aggregate reproductive outcomes and sex-specific fertility liabilities from genotypic information.

Heavy metals act as essential metalloprotein cofactors in numerous physiological processes but can become toxic when non-essential metals accumulate or when essential metals are in excess. As plants continuously encounter heavy metals through their roots, they have evolved complex homeostatic mechanisms to regulate metal uptake and distribution. The Heavy Metal ATPase (HMA) gene family encodes a group of heavy metal transporting P-type ATPases that have been linked to stress resistance and nutrient supply. Here, we used a bioinformatics approach to identify and characterise 13 HMA genes containing characteristic P1B-type ATPase domains and motifs in the barley Morex V3 reference genome. The genes are located on five of the seven barley chromosomes. Phylogenetic analysis revealed that they cluster into five sub-clades, including one clade unique to barley. Expression profiling across multiple datasets showed distinct temporal and tissue-specific expression patterns among HvHMAs, with several members exhibiting significant transcriptional responses to specific biotic and abiotic stresses. By utilising recently available pan-transcriptomic and pan-genomic resources, we have identified substantial allelic diversity and inter-accession variation in HvHMAs. Our findings suggest that HvHMAs have functions extending beyond canonical heavy metal homeostasis and warrant further investigation for their potential roles in broader physiological and stress-related processes.

Drosophila suzukii (Matsumura, 1931, Diptera: Drosophilidae) is a globally invasive pest of soft-skinned fruits that is currently controlled largely through the use of broad-spectrum insecticides. Increasing resistance to pesticides and regulatory pressures have motivated the development of genetic control strategies. We previously developed a CRISPR/Cas9-based homing gene drive targeting the coding sequence of the female-specific exon of the sex-determination gene doublesex, achieving highly efficient inheritance (94-99%) in both male and female germlines. A major limitation of homing gene drives is the formation of resistant alleles that evade cleavage yet retain gene function. Multiplexing guide RNAs (gRNAs) could reduce the formation of such functional resistance alleles. Here, we generated and tested homing constructs expressing one, two, or three gRNAs targeting different regions of the female-specific exon of doublesex, including a splice-junction target site. A single gRNA targeting the splice junction supported high inheritance in males but showed reduced efficiency in females. Combining this gRNA with a coding sequence-targeting guide further reduced drive efficiency, particularly in the female germline. Constructs expressing two gRNAs performed similarly whether guides were linked by transfer RNA (tRNA) sequences or expressed from independent promoters. Constructs expressing three gRNAs using tRNA processing showed consistently low drive inheritance in both sexes, likely a consequence of reduced cleavage efficiency due to inefficient gRNA production. Inheritance was significantly higher in male than female germlines for several constructs, indicating that germline context strongly influences drive performance. Our findings highlight that the approach used for multi-gRNA expression, target site choice and sex-specific germline environments can impact gene drive efficiency, and emphasize the need to optimize construct design within the biological context of the target species.

Entamoeba histolytica is a clinically important pathogenic eukaryote and the causative agent of amoebic dysentery. Entamoeba dispar, a nonpathogenic commensal species that resides in the human colon, is the closest sibling species, and serves as an appropriate comparator for genome-wide analysis. Although the genome of E. histolytica is approximately 26.9 Mb, and the largest known genome within the genus, that of E. invadens, is approximately 40.9 Mb, obtaining high-quality assemblies in this genus has remained challenging due to extensive repetitive regions, tRNA gene arrays, and aneuploidy. Here, we used PacBio HiFi sequencing to assemble the genomes of the pathogenic E. histolytica and the nonpathogenic E. dispar. We reconstructed all 36 chromosomes of E. histolytica and 35 chromosomes of E. dispar, assembling each as a single continuous DNA sequence (contig). The two species exhibited high genome-wide nucleotide similarity and conserved synteny at the amino acid level. At one end of each chromosome, we identified tRNA arrays, whereas the opposite end lacked such arrays, resulting in an asymmetric chromosomal architecture. Analysis of unique-read depth revealed widespread aneuploidy in both species: E. histolytica is predominantly tetraploid, whereas E. dispar is diploid, a conclusion further supported by SNP allele-frequency distributions. These assemblies provide a robust foundation for comparative genomics in Entamoeba and offer detailed insights into chromosome-end structure and ploidy.

Rice (Oryza sativa) is highly sensitive to salinity, yet the metabolic mechanisms underlying salt tolerance remains incompletely understood. In this study, we performed leaf tissue-specific untargeted metabolomic profiling of the salt-tolerant introgression line JN100 (JN), its donor parent Nona Bokra (NB), and its recurrent parent Jupiter (JU) to characterize metabolic responses to salt stress. Comparative analysis identified differentially accumulated metabolites (DAMs) spanning diverse chemical classes, including amino acids, sugars and carbohydrates, lipids, organic acids, cofactors, electron carriers, and nucleotides. Under salt stress (SS), 201 DAMs (89 upregulated and 112 downregulated) were detected in JN relative to JU. Notably, metabolites such as allantoin, glycitin, nicotinamide ribotide, D-arabinono-1,4-lactone, violanthin, L-methionine S-oxide, ribitol, lysine, rutin, glutamine, pantothenic acid, and quinic acid, showed significant differential accumulation. Pathway enrichment analysis revealed significant enrichment of arginine biosynthesis, purine metabolism, and alanine, aspartate, and glutamate metabolism, indicating extensive reprogramming of nitrogen and energy-associated metabolic pathways under salinity stress. Integration of transcriptomic and metabolomic datasets from the SS experiments further identified ten differentially expressed genes (DEGs) associated with the metabolite network in the JN vs. JU comparison. Among these, OsDHQDT/SDH, OsFd-GOGAT, phenylalanyl-tRNA synthetase, OsP5CS1, OsP5CS2, and a pyridoxal phosphate-dependent transferase were linked to metabolites involved in shikimate, amino acid, and proline metabolism. Collectively, these results demonstrate that salinity tolerance in rice is associated with coordinated transcriptional and metabolic reprogramming that supports oxidative stress mitigation and adaptive stress responses.

In plant breeding and genetics, recent advances in high-throughput phenotyping are beginning to meet the growing demand for large-scale, high-quality phenotypic data that emerged after the development of next-generation sequencing technologies. Recent developments in phenomics have been incorporated into almond breeding programs, facilitating the large-scale acquisition of quantitative phenotypes and the dissection of the genetic architecture underlying morphological and quality-related traits. The implementation of a high-throughput phenotyping platform integrating RGB and hyperspectral imaging with genotyping using the 60K almond SNP array enabled the large-scale characterization of almond populations and the identification of 567 robust marker-trait associations across 66 traits. These analyses revealed two major genomic hotspots on chromosomes 2 and 5 associated with morphological and quality-related traits. These regions harbored biologically relevant candidate genes, including genes associated with OVATE family proteins, brassinosteroid signaling, protein ubiquitination, and acyl-CoA metabolism, as well as other regulators of organ growth, cell proliferation, hormone signaling, and seed development. Furthermore, a novel candidate gene encoding a COMT-like O-methyltransferase involved in lignin biosynthesis was identified and proposed to contribute to shell hardness, a major genetically controlled trait in almond. Together, these findings demonstrate the potential of integrating high-throughput phenomics and genomics to dissect complex traits, identify candidate genes, and accelerate genomics-informed breeding in almond.

Animals fight for resources to obtain fitness benefits; most contests are intrasexual, and males tend to fight more than females. Although the genetic basis of male aggression is well studied, we know little about the genetic variation of female aggression. Female aggression varies with reproductive status and is potentially influenced not only by her genotype, but also by the genotype of her mate. Here we measured both male and female aggression in a set of Drosophila melanogaster inbred lines by competing each line against a standard competitor. Aggression varied among lines for both sexes, but male and female aggression were not correlated. Female aggression for many lines increased with mating, as expected, but not all lines changed aggression. However, when females were mated to males of different lines, male genotype did not affect the post-mating change in aggression, suggesting that ejaculate-mediated effects do not vary across these lines. The aggression level of the standard opponent was positively correlated with that of focal individuals indicating that individuals modulate their behavior according to the genotype of their opponent.

The potato leafhopper, Empoasca fabae (Harris, 1841), is a highly polyphagous, migratory insect pest of eastern North America that feeds on more than 200 herbaceous and woody plant species, causing substantial losses to forage and field crops. Despite its agricultural and ecological importance, no genome has been available for this species. Here, we present the first chromosome-level genome assembly of E. fabae, generated from Oxford Nanopore long reads, Illumina short reads, and Omni-C proximity-ligation data. The final assembly spans 908 Mb across 132 scaffolds, with 99.8% of the assembly captured in ten chromosome-length scaffolds (nine autosomes and an X chromosome) with a scaffold N50 of 96.2 Mb. The assembly is highly complete, recovering 92.4% of conserved hemipteran single-copy orthologs, and is composed of 47.6% repetitive sequence, dominated by long terminal repeat retrotransposons and unclassified elements. Read-depth comparison between male and female individuals supports assignment of a single sex-linked chromosome, consistent with an XO sex-determination system. BRAKER3 gene annotation predicted 31,406 protein-coding genes after retaining the longest isoform per locus. Comparative genome analysis against the two closest related Typhlocybinae species with genomes available, Matsumurasca onukii and Hebata decipiens, revealed extensive chromosome-scale collinearity, while defining a shared core gene repertoire. This reference genome provides a foundation for comparative and population genomic studies and for investigating genetic traits in this economically important crop pest species.

Celiac disease (CeD) is an immune-mediated condition triggered by wheat gluten in genetically predisposed individuals. The immune reaction in people with CeD is driven by particular gluten amino acid sequences, or immunogenic epitopes. Some of these epitopes elicit strong immune responses in the majority of CeD patients and are designated as immunodominant epitopes. Previous research has shown correlations between the amount of immunogenic wheat epitopes consumed and the onset of CeD, suggesting that reducing wheat immunogenic epitopes may reduce CeD incidence at the population level. Gluten consists of gliadins and glutenins, with gliadins having the majority of the immunodominant epitopes and glutenins playing a major role in dough strength and breadmaking quality (BMQ). This study used radiation-induced deletions, chemical mutagenesis, and natural variation in wheat (Triticum aestivum) to generate genetic stocks with reduced immunogenic epitope content. Most lines were developed in the wheat cultivar Summit, for which we produced a full genome assembly and annotation. We used exome capture to characterize these deletions and identify prolamins located within and outside the deletions. We combined different deletions and developed molecular markers to facilitate their deployment. For chromosome arms 1BS and 1DS, we generated two alternative lines: one lacking immunogenic epitopes for the development of CeD-safe genetic stocks for research purposes, and another retaining selected glutenins for breeding commercial lines with reduced immunogenicity and adequate BMQ. By making these non-transgenic genetic stocks publicly available, we aim to accelerate the development of wheat varieties with reduced immunogenicity and, eventually, a fully CeD-safe wheat.

The availability of pangenome and resequencing of wheat collections have facilitated the discovery of gene-trait associations in wheat. Yellow stripe-like (YSL) proteins play a key role in the uptake and translocation of metals and yet have not been fully identified and analyzed at the genome-wide level in wheat. In this study, 26 TaYSL genes were identified and divided into four distinct clades, each clade sharing similar domains and motif compositions. Most genes were upregulated under iron deficiency, whereas homoeologs of TaYSL1 were downregulated. Both SNP-based and haplotype-based association studies were used to dissect the role of TaYSLs underpinning grain iron contents (GFeC) and zinc contents (GZnC) in wheat. TaYSL6-2B and TaYSL16-1A haplotypes showed strong association with GFeC, and TaYSL14-6A showed strong association with GZnC in multiple field trials. The distribution of favorable haplotypes in global wheat collection of [~]3000 accessions showed that majority of haplotypes were more prevalent in landraces and winter wheat compared to modern cultivars and spring types, indicating their potential for use in breeding. The combination of favorable haplotypes of three YSL genes associated with GFeC and GZnC were very rare, and most of the wheat accessions has single or double favorable haplotypes. These findings provide the first comprehensive characterization of the TaYSL gene family in wheat and identify significant SNPs and elite haplotypes that can be utilized for genetic improvement and biofortification.

The vine mealybug, Planococcus ficus, is a globally invasive pest of grapevine and a vector of leafroll viruses. Like other mealybugs, it reproduces through paternal genome elimination, a sex-determination system that operates without sex chromosomes and is associated with extreme sexual dimorphism. To characterize genome organization and sex-biased expression in this species, we generated a long-read reference genome spanning 369 Mb with 23,489 annotated genes and macrosynteny conserved with the citrus mealybug, Planococcus citri. Resequencing of four California field individuals yielded a first whole-genome estimate of nucleotide diversity and 132 microsatellite markers for population monitoring. Among 2,129 candidate secreted proteins, a conserved core is shared with P. citri, but each species carries a distinct set of lineage-specific effectors. Comparing adult male and female transcriptomes, we found sex-biased expression to be pervasive and skewed toward females: 41% of tested genes differed between the sexes, with female-biased genes both more numerous and showing larger fold changes. These female-biased genes were not randomly distributed but concentrated in discrete blocks of coordinately expressed, tandemly duplicated gene families, a pattern not previously described in a mealybug. Male- and female-biased secreted proteins also differed in origin, with male-biased proteins drawn from a conserved repertoire shared with P. citri and female-biased proteins spanning a more lineage-specific pool. Together, these results reveal a female-skewed, spatially clustered architecture of sex-biased expression in a mealybug that lacks sex chromosomes, and provide genomic resources for managing an invasive vineyard pest.

Accessory genes are thought to contribute to fungal adaptation and pathogenicity by modulating host immunity, while core genes play crucial roles in maintaining fundamental biological processes. Rust fungi (order Pucciniales) are obligate biotrophic plant-pathogens and infect economically relevant crops. Here, we characterize core and accessory gene repertoires across rust fungi, with a particular focus on Phakopsora pachyrhizi, the causal agent of Asian soybean rust. Across Pucciniales genomes, accessory genes represented the largest fraction of gene content (~44.6% on average), whereas core genes accounted for a smaller proportion (~18-35%). Notably, variations in accessory gene content among rust fungi are perhaps attributed to lineage-specific gene expansions and losses. Core gene content was positively correlated with total gene number across Pucciniales genomes, suggesting retention after gene duplication events, consistent with their essential biological functions. Among P. pachyrhizi genes expressed during soybean infection, core effectors were associated with cysteine-rich proteins, pectin-degrading enzymes, and SPFH/Band 7 family, while accessory effectors included phosphatidylethanolamine-binding proteins, trehalose phosphatases, and CFEM domain-containing proteins. The in-plant induced core and accessory genes in P. pachyrhizi also comprised multiple families of CAZymes (GH5/GH7 cellulases, CE5 cutinases, CE8 pectinesterases, CE4/GH18 chitin-modifying enzymes); proteases (aspartyl proteases, serine carboxypeptidases, alpha/beta hydrolases); transporters (amino acid permeases, ferric reductase-like transmembrane proteins, and OPT oligopeptide transporter), and transcription factors (bZIP, GATA zinc finger, STE-like, and homeobox KN). Our study highlights that core and accessory gene families have shaped P. pachyrhizi-soybean interactions, identifying promising targets for functional studies aimed at elucidating host-adaptation mechanisms in rust fungi.

Potato leafroll virus (PLRV) is an economically important viral disease of potato (S. tuberosum). Genetic resistance to this phloem-limited virus is rare, and no cloned resistance (R) genes have been reported. Rladg confers resistance to PLRV in an Andean potato landrace, LOP-868 (Velasquez et al. 2007). We identified the functional Rladg gene as a homolog of the tomato TIR-NLR-encoding Bs4. Rladg interacts with the serine protease domain of the PLRV protein P1, which is essential for virus replication. This recognition is independent of the proteases enzymatic activity, and the Rladg immune receptor oligomerizes upon direct association with the protease. Like PLRV, many poleroviruses contain a serine protease. Despite their diverse amino acid sequences, these proteases are predicted to share similar structures. Rladg recognizes all ten tested polerovirus proteases, suggesting a conserved structural recognition mechanism. We propose that Rladgs broad recognition capacity could enable resistance to poleroviruses in many crop species. Rladg is the first R-gene reported to confer resistance to a phloem-limited pathogen and could provide enhanced resistance to many economically important poleroviruses.